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We demonstrate on examples that a simultaneous study of the Bose- 
Einstein correlation function and the invariant momentum distribution 
can be very useful in distinguishing various hydrodynamic models, 
which describe separately the short-range correlations in high energy 
hadronic reactions as measured by the NA22 collaboration. We also 
analyze Bose-Einstein correlation functions, measured by the NA44 
experiment at CERN SPS, in the context of the core-halo model. Val- 
ues for the core radius and the fraction of direct bosons are obtained, 
and found to be independent of the structure of the correlation function 
at small relative momenta of Q < 40 MeV. 



1 Introduction 

The NA22 Collaboration performed recently a detailed study of Bose-Einstein Cor- 
relation Functions (BECF-s) in one- two- and three-dimension£l, for tt + / K + + p re- 
actions at 250 GeV. Their study concluded as follows: "Our data do not confirm 
the expectation from the string- type model ... . A good description of our data is 
... achieved in the framework of the hydrodynamic expanding source modelfli'ii'iiH 
... . Alternatively, our data are also described in a non-expanding, surface-emitting 
fireball-like sources^... ." The above cited hydrodynamic models, successful in 
describing the NA22 data, belong to the following three different classes: i) Non- 
expanding, spherically symmetric fireball ii) Expanding, spherically symmetric 
shellfll. Hi) Longitudinally expanding, cylindrically symmetric models with possible 
transverse flow, transverse and temporal temperature profileJEi'ii. 



2 Particle correlations and spectra in various hydrodynamic models 

The experimental evidence reported in ref. I indicates that it is rather difficult to 
distinguish the Fourier-transformed emission functions of the models of type i) - in). 
However, the same emission function, which determines the two-particle BECF-s, 
prescribes also the single-particle spectra. In the forthcoming, we shall summarize 
the particle spectra and the BECF-s of the models and discuss what are the essential 



differences among these. We shall utilize the Wigner-function formalism along the 
lines of refs@0i'i, where the particle emission is described by the emission function 



S(x,p), x = (t, r) and p = (E, p) with E = y |p| 2 + rn 2 . The two-particle BECF for 
completely chaotic sources is givei0lll0'i by 



(n(n - 1)> JVi( Pl ) Af!(p 2 ) |S(0,A')| 2 ' 

where N(p) = EdN/dp = S(0,p) stands for the single-particle inclusive invariant 
momentum distribution, normalized to the mean multiplicity as / N(p) dp/ E = (n) 
and S(Ak, K) = J d A x S(x, K) exp(iA/c • x) with Ak = p% — p 2 , and K — {p\ + p 2 )/2, 
where Ak ■ x stands for the inner product of the four-vectors. 

2. 1 Spectra and Correlations for Kopylov and Podgoretskii model 

The source considered by Kopylov and Podgoretskii (KP) in refslE was one of the 
earliest models of particle emission in multiparticle physics: a uniformly illuminated, 
surface emitting sphere. The emission function itself was not discussed in detail, since 
the source was assumed to be non- expanding, static. If we assume that the source is 
thermalized, the KP emission function is given by 

S KP (x,p) = j^- s E c S(Rkp - H) eM-t/TKP) eME ' /T) _ r (2) 

which is defined in the CMS of the source. The degeneracy factor is indicated by g, 
the pre-factor E c is due to the invariant normalization of the momentum distribution. 
The parameter Rkp is the radius of the source, r^p is the decay time of the quanta, 
or in another interpretation, the thickness of the emitting layer, E c is the energy of 
the particles in the CMS of the source and T is the temperature of the fireball. The 
single-particle spectrum can be obtained as 

dN KP , N g 1 

E-t- = N KP (p) = -^E C V KP — — 3 

dp (27m) exp(£/ c /i ) — 1 

where Vkp — ^R% p tkp is the effective volume of the fireball. In Boltzmann ap- 
proximation, this invariant momentum distribution can be rewritten as 



dN g f m t \ 

N Kp{p) = 1 — j- ^ jz ., 3 mt cosh(y - y )V KP exp -- — — , (4) 

2irm t dm t dy (27ih) \ T KP [y) ' 

T K p(y) = T/cosh(y -y ), (5) 



where m t and m t = ^Jp 2 . + p 2 + m 2 is the transverse mass, y = 0.5 \og((E + p z ) / (E — 
p z )) denotes the rapidity of the particle, and y stands for the rapidity of the CMS of 
the fireball. The effective rapidity-dependent temperature distribution , as given by 



eq. decreases fast with increasing difference between the rapidity of the particle 
as compared to the CMS of the fireball. This 1/ cosh(y — y ) decrease is a typical 
result for non-expanding thermalized fireballs, as discussed in ref@. The spectrum of 
eq. (0) can be re-written approximately as 

dN g ( (y-y ) 2 \ 

2nm t dm t dy " j2wKf m t cosh(y - y Q )V K p exp ^ iKrftT ) ' (6) 

where Arjj, = T/m t . This yields a specific transverse mass dependence for the 
rapidity- width of the spectrum at a given value of m t , which can be checked easily in 
experimental data analysis. The BECF for the KP emission function is given as 

C(q T ,K) = l + \[I(R KP q T )] 2 (l + r 2 KP q 2 )-\ (7) 

where I(x) = 2J 1 (x)/x and J±(x) is the Bessel function of first order, with 1(0) = 1 
and I(oo) = 0. The parameter A is introduced phenomenologically to account for 
intercepts 1 + A < 2, = Ak • K/|K| is the relative momentum transverse to K and 
go = E\ — E2 is the energy difference in the CMS. This BECF is independent of the 
mean momentum K due to the non-expanding source described by the KP model. 



2.2 Spectra and Correlations for Expanding Shells 

The model introduced by S. Pratt (P) and discussed recently by J. Bjorken corre- 
sponds to a uniformly illuminated, expanding shell! given by 

S p( x 'P) = 777^3 E C 5(R P - \r\)exp(-t 2 /r 2 P )- ' 



(27rhy ' ^ exp(p-u(|r|)/T)-l' 

which is defined in the CMS of the source. The parameter Rp is the radius of the 
source, Tp is the decay time of the quanta or the thickness of the emitting layer 
in the CMS of the source and w(r) = 7(1, vr/Rp) is a spherically symmetric flow 
profile characterizing the expansion of the shell, with 7 = 1/y/l — v 2 . This model 
corresponds to a three-dimensional, spherically symmetric expansion in the rest frame 
of the source, and the KP model is recovered in the v — limiting case. 

The single-particle spectrum can be evaluated in Boltzmann approximation as 

dN g sinh(a) 
~ -3 m t cosh(y - y )V P exp(-m t /T P (y)) (9) 



2'Km t dm t dy (2irh) ' a 

T P (y) = -[1 + CV7 2 )] (10) 

7 cosh(y - y ) 

where Vp = AnRp^Trrp) 2 is the volume of the fireball, a = a(p) = vj\p\/T is a 
momentum-dependent parameter related to v , the surface velocity of the shell and 
Tp(y) is the effective rapidity - dependent slope-parameter distribution. Thus for 



small values of the parameter a(p) the effective temperature is decreased as com- 
pared to the static fireball (KP) case, however, the rapidity-width of this effective 
temperature distribution remains unchanged even if a small spherical expansion is 
included. The BECF for the spherically expanding shells can be written in a direc- 
tional dependent, analytic form as given by Eq. (13) of refi. Although this analytic 
form is too complicated to be repeated here, one arrives at a simpler expression^, 

(C(q,K)) = l + Aexp(-0.5i% p(? 2 ) exp(-0.5 r 2 g 2 ), (11) 

after averaging over the direction of the relative momentum in the CMS of the source. 
Here q = |Ak| in the CMS and the radius parameter read@ as 

R %P = J R P [(a(p)tanh(a(p))^ 1 -sinh(a(p))- 2 ] 1/2 - (12) 

The momentum dependence of the effective radius parameter R*,p is a direct conse- 
quence of the expansion of the source. 



2. 3 Spectra and Correlations for Longitudinally Expanding Systems 

Systems which are dominantly expanding longitudinally appear to be the relevant 
models to high energy heavy ion reactions at CERN SPS energy region, y/s = 
20GeVA. The detailed presentation and elaboration of these typ e of models is outside 
of the scope of the present contribution, we recommend refsilBii for further details. 
However, we summarize here approximate results for the single-particle spectra: 



dN 



■ m t cosh(y - 77 s (p))V*(p) exp(-m t /T L (y)), (13) 



27im t dm t dy (2nh) 3 

2nm d t dm t dy " (2^™*°°^ ~ ^{~^^) ' (14) 

K(p) ^ V (T/m t ) a -\ (15) 

ar(y-yof 



T l(v) = T—n 1 . . x2 » and A 2 y(m t ) = A 2 r] + T/m t . (16) 



These relations indicate that the rapidity dependence of the temperature parameter 
Ti(y) can be described with a new fit parameter a? (which was showri to be related 
to Arj, the total longitudinal extension of the particle emitting source, e.g. = 
for infinite systems). A more direct access to the longitudinal size of this expanding 
system is provided by the transverse mass dependence of the rapidity-width of the 
invariant momentum distributiori'i, because A 2 y(m t ) ~ A 2 rj for mt >> T. 
The BECF-s for such systems can be written to the following form: 

C(Ak, K) = 1 + A*(K) exp (-R 2 slde {K)Q 2 slde - R 2 out mQ 2 out - R 2 ong (K)Q 2 ong ) x 

exp (-R 2 olong (K)Q out Q long ) , (17) 



where the intercept parameter can be interprete in the core-halo model, the side, 
out and longitudinal radius parameters as well as the out-long cross-terni may in 
general depend on the mean momentum K in a complicated manneriS!. However, in 
some specific limiting case J this can be simplified as R s ide — Rout — Rimg oc 1/ yjrh~t . 

2.4 How to distinguish the three model-classes? 

The following tests can be performed experimentally to check the hypothesis of Kopy- 
lov and Podgoretskii, Pratt and Bjorken or the model-class of longitudinally expand- 
ing, finite systems in a more detailed manner: 

1. Measure the m t dependence of the effective rapidity-width of the N(p) dis- 
tribution at a fixed value of m t , and try to fit the result with Eq. (|16|). For 
static or spherically expanding shells, the parameters Ay(m t ) decrease to as 
l/ A /mt, while for longitudinally expanding finite systems the large mt limit is 
Ay(oo) = At] > 0, the longitudinal size of the system in space-time rapidity. 

2. Measure the rapidity dependence of the temperature parameter of the N(p) 
distribution at a fixed value of y, and try to fit the result with Eq. fllEp . A 
slow drop of the effective temperature with increasing values of \y — yo\ results 
in the breakdown of the static fireball picture. Longitudinally very extended, 
expanding systems are predicted to have a Lorentzian effective temperature 
distribution, T e ff(y) = T*/(l + axiy — yo) 2 ) where the parameter clt carries 
information about the longitudinal extension of the sourcJl 

3. Check if the parameters Rkp, Rp,Tp and r^p of the Bose-Einstein correlation 
function of eq. (|7|) are independent or not of various values of K. Static models 
predict Rkp(K) = const and r^p(K) = const, a deviation from this behav- 
ior results in the breakdown of the KP or any other static model. Spherically 
symmetric expanding shells also have a specific transverse momentum depen- 
dent R*^p(K) parameter, decreasing with increasing values of |K| as given by 
eq. (0). Longitudinally expanding finite systems have momentum dependent 
longitudinal radius parameter, the transverse radius parameters may or may 
not be transverse mass dependent. 

Let us now focus on the interpretation of the intercept parameter of the BECF-s. 
3 Core-halo model analysis of NA44 data 

Having first performed a Monte Carlo simulation to justify the analysis technique to 
be used, we then analyze the Bose-Einstein correlation functions from NA44^1'lil. 



3.1 Monte-Carlo test of the method 



In a Monte-Carlo study, we simulate the actual and background distribution (A(Q) 
and B(Q)) of particle pairs. For the CERN experiment NA44, studying S+Pb colli- 
sions, the shape of B(Q) is given approximately by B(Q) = Q 3 e~( 3 ' 6 ® >, which form 
reproduces the experimentally measured background distribution^. This background 
distribution is peaked at 25 MeV similarly to the NA44 data. We have started with 
a core-halo model correlation functional, assuming Gaussian source functions for the 
core and for the halo, choosing the core radius parameter to be R c = 4 fm. The 
radius parameter of the halo was taken to be Rh = 40 fm and the core fraction was 
f c = 0.7. The actual distribution A(Q) was sampled according to the formula 

A(Q) = B(Q)[l + f^exp(-RlQ 2 ) + 

2 f c (1 - f c ) exp(-0.5(i? c 2 + R 2 h )Q 2 ) + (1 - f c f exp(-R 2 h Q 2 )] (18) 

which corresponds to the full correlation function including correlations of (c, c), (c, h) 
and (h, h) type of particle pairs. Here c refers to the core, which is assumed to be 
resolvable since R c < fi/Qbi n , where the two-particle relative momentum resolution 
Qbin is about the size of a bin (cca. 10 MeV for NA44). In the core-halo mode@, the 
halo (index h) is assumed to change on large length-scales, which are un-resolvable 
by the two-particle correlation measurements, Rh > h/Q min . 

Having sampled the actual and background distributions of particle pairs in a 
Monte-Carlo simulation, the correlation function is now calculated in each bin of size 
Qbin as C 2 (Q) = b(o\- This form is fitted by the expression 

C 2 (Q) = 1 + A* exp(-RlQ 2 ). (19) 

If the core-halo model is applicable, it predicts that A* = f 2 and i?* = R c after the 
fitting, i.e. the intercept can be utilized to measure the fraction of particles emitted 
by the core, and R* the radius parameter will coincide with the radius parameter 
characterizing the core. Furthermore, if the assumptions of the core-halo model are 
satisfied by some data, then the fitted A* and i?* parameters should become indepen- 
dent of the fitted region, obtained by excluding more and more part of the correlation 
function at the lowest values of Q, as long as the inequalities Ti/R h < Q m i n < Ti/R c 
are satisfied by Qmin- We chose the range from Qu n = 10 < Q m %n < 60 MeV. 

We have checked in the Monte-Carlo simulation if the effect of the halo on the 
deviation from a Gaussian form does really cancel or not from the simulation by 
increasing the value of Q m in, the size of the low Q region excluded from the fit. We 
have found that indeed A* = f 2 ~ 0.5 and R* = R c = 4.0 fm was reproduced by 
the fitting within the statistical uncertainties even for as large values of the excluded 
region as Q m in = 60 MeV, when already more then one half of the peak was removed 
from the correlation function, utilizing 300 K pairs. Armed with this experience, we 
then applied the method to the analysis of NA44 data for pions and kaons in various 
one-dimensional slices of the correlation functions at given values of K. 



3.2 Analysis of NA44 data 



In order to extract values for R andf c from data, we first fitted the side, out and 
longitudinal slices of the NA44 datai^l on S + Pb reactions at 200 AGeV for pions 
at low and high transverse mass as well as for kaons. Thus we obtained radius 
parameters and intercept parameters (with their errors, respectively), as a function 
of Qmin- We determined these fit parameters in the Q m in = — 60 MeV region. The 
fitted R*(Qmin) and X*(Q m in) parameters up to Q m in = 40 MeV/c were (meta)-fitted 
with a Qmin independent constant. Our findings are summarized in Table 1. More 
details of the study shall be reported elsewheriB. 



parameter 


high p t 7T + 


low Pt 7T + 


K + 


R c ou t (fm) 

Rc side (fm) 
Rc long (fm) 
fc out 
fc side 
fc long 


2.92 ± .13 
2.90 ± .18 
3.31 ± .16 
.704 ± .012 
.735 ± .021 
.738 ± .014 


4.29 ± .13 
4.24 ± .26 
5.43 ± .30 
.725 ± .011 
.647 ± .021 
.724 ± .015 


2.54 ± .18 
2.22 ± .19 
2.67 ± .22 
.802 ± .027 
.736 ± .033 
.789 ± .029 



Table 1: Extracted values for R c and f c from NA44 S+Pb data. 

In ref@ the halo has been interpreted as being created essentially by the decay 
products of the u, r\ and r( resonances. The measured core fractions of Table 1 are 
indeed in the vicinity of the core fractions of Fritiof and RQME0 if cu, f] and 7/ are 
taken as unresolved long lived resonances, as indicated by Table 0. Thus we do not 



fc,NA44 (high pt 7T+ ) 


fc,NA44 (l0W p t 7T+ ) 


fc, Fritiof 


fc,RQMD 


0.71 ± 0.01 


0.71 ± 0.01 


0.75 


0.80 



Table 2: Core fractions f c from NA44 S+Pb data compared to core fractions from RQMD 
and Fritiof, when u>, n and 7/ are taken as unresolved long-lived resonances 17 . The core 
fractions for the NA44 low pt and high pt sample are obtained from a simultaneous fit to 
(fc,out, fc,side, fc,iong) by a constant value, and the result is rounded to two decimal digits. 

find a complicated, non-Gaussian structure of the measured BECF in S+Pb reactions 
at CERN SPS, in contrast to some earlier expectations which argued that resonance 
decays, coupled to a freeze-out hyper-surface obtained from (hydro) dynamic evolution 
may lead to resolvable deviations from a Gaussian structure. Such a deviation has 
been predicted by detailed simulations of resonance decay effects by calculations with 
SPACEF0, HYLANDEB0, and recently by H. Heiselber|l. Assuming that such 
non- Gaussian structures are created by the effect of r],r]' and uj resonance decays to 
pions, but that they are experimentally not yet resolvable in the NA44 S + Pb — > 
2tt + X experiment, we obtain a consistent interpretation of tl^data and a physical 
interpretation of the measured intercept and radius parameterfi2Hi-lu. 
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